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The host galaxy of a narrow-hne Seyfert 1 galaxy RE 
J1034+396 with X-ray quasi-periodic oscillations 
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ABSTRACT 

Using simple stellar population synthesis, we model the bulge stellar contribution in 
the optieal spectrum of a narrow-line Seyfert 1 galaxy RE J1034-I-396. We find that 
its bulge stellar velocity dispersion is 67.7 ± 8 km . The supermassivc black hole 
(SMBH) mass is about (1 — 4) x IO^Mq if it follows the well-known Mbh — ct* relation 
found in quiescent galaxies. We also derive the SMBH mass from the H/3 second 
moment, which is consistent with that from its bulge stellar velocity dispersion. The 
SMBH mass of (1 — 4) x IO^Mq implies that the X-ray quasi-periodic oscillation 
(QPO) of RE J1034-F396 can be scaled to a high-frequency QPO at 27-108 Hz found 
in Galactic black hole binaries with a 10 Mq black hole. With the mass distribution 
in different age stellar populations, we find that the mean specific star formation 
rate (SSFR) over past 0.1 Gyr is 0.0163 ± 0.0011 Gyr'^ the stellar mass in the 
logarithm is 10.155 ± 0.06 in units of solar mass, and the current star formation rate 
is 0.23 ± 0.016 Mq yr"!. RE J1034+396 does not follow the relation between the 
Eddington ratio and the SSFR suggested by Chen et al., although a larger scatter in 
their relation. We also suggest that about 7.0% of the total Ha luminosity and 50% 
of the total [O ii] luminosity come from the star formation process. 

Key words: galaxies: bulge — galaxies: nuclei — black hole physics — galaxies: 
stellar content — Galaxy:individual: RE J1034-I-396 
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1 INTRODUCTION 

It is thought that active galactic nuclei (AGN) are scaled-up 
versions of Galactic black hole binaries (BHBs; Gierlinski et 
al. 2008 and reference therein). The power spectrum of the 
X-ray variability in BHBs shows the quasi-periodic oscilla- 
tions of 0.01-450 Hz(QPOs; Rimillard & McClintock 2006). 
It is believed that we can also find the QPOs in AGN, on 
the similar behavior of accretion flow around the black hole 
(BH). Comparing to a BH of 10 Mq in BHBs, the typi- 
cal high-frequency QPOs of 100 Hz would be smaller by 
10® (i.e., ~ 10~* Hz) for a supermassive black hole (Mbh 
; SMBH) of IO^Mq if the QPO frequencies scale inversely 
with the BH mass (Rimillard & McClintock 2006). 

With a long XMM-Newton observation (91 ks), Gier- 
linski et al. (2008) detected a significant QPO signal [v = 
2.7 X 10~*Hz, corresponding to a period of about 1 hour) 
for a nearby (z=0.043) spiral active galaxy RE J1034-I-396 
(J2000, RA=158.66082, DEC=39. 64119). It is optically clas- 
sified as a narrow-line Seyfert 1 galaxy (NLSl), with similar 
small line width for the high-ionization and low-ionization 
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emission lines (Puchnarewicz et al. 1998), and a very soft 
X-ray spectrum (Grupe et al. 2004; Casebeer et al. 2006). 
NLSls are thought to be a special subclass of AGN har- 
boring relatively small but growing SMBHs, compared with 
other broad-line Seyfert 1 galaxies (BLSls; e.g. Osterbrock 
& Pogge 1985; Boiler, Brandt & Fink 1996; Mathur 2000). 
By methods of the standard accretion disk model fitting 
of spectral energy distribution (SED; Puchnarewicz et al. 
2001), the slim disk fitting of SED (Wang et al., 1999; Wang 
& Netzer 2003), the full width at half-maximum (FWHM) 
of the H/3 line, [O ill] FWHM, and soft X-ray luminosity, 
the SMBH mass in RE J 1034-1-396 can not be determined 
very weU, from 6.3 x 10^ to 3.6 x US' Mq (Bian et al. 2004; 
Wang & Lu 2001). Its QPO type (high-frequency or low- 
frequency QPOs) can not be uniquely identified (Gierlinski 
et al. 2008). 

In the past two decades, there has been striking progress 
in finding more reliable methods to calculate SMBHs masses 
in AGN through the line width, AU, of H/? (or Hq, Mg ll , 
C IV ) from the broad line region (BLR) and the BLR size, 
-Rblr (e.g., Kaspi et al. 2000; McLure & Dunlop 2004; Bian 
& Zhao 2004; Peterson et al. 2004; Greene & Ho 2005b; Bian 
et al. 2008). There are mainly two ways to parameterize the 
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line widths of broad emission lines, i.e., FWHM and the 
second moment (crjijie) (e.g. Bian et al. 2008). The second 
moment of the broad components of the H/3 line provides a 
more precise measurement of the SMBH mass because the 
H/3 profile is non-Gaussian in NLSls. It is suggested that 
the mean value of the SMBH masses in NLSls from the 
H/3 second moment is larger by about 0.50 dex than that 
from FWHM (Bian et al. 2008). The well-known Mbh - cr, 
relation of inactive galaxies (Tremaine et al. 2002) can also 
be used to estimate the SMBH mass when the stellar velocity 
dispersion, a,, is available (Kauffmann et al. 2003; Onken 
et al. 2004; Greene & Ho 2006). 

The X-ray QPOs frequency may depend on the mass 
and the spin of BH (Rimillard & McClintock 2006). The 
SMBH masses are important in the study of QPOs in AGN. 
Here we use the the Sloan Digital Sky Survey (SDSS; Abaza- 
jian et al. 2009) data of this object to determine its SMBH 
mass by its cr, and H/3 second moment instead of H/3 FWHM 
in Bian et al. (2004), as well as the star formation history of 
its host bulge by the simple stellar population (SSP) synthe- 
sis. All of the cosmological calculations in this paper assume 
Ho = 70 km s~^ Mpc~\ = 0.3, and = 0.7. 



2 DATA AND ANALYSIS 



The SDSS used a dedicated 2.5-m wide-field telescope at 
Apache Point Observatory near Sacramento Peak in South- 
ern New Mexico to conduct an imaging and spectroscopic 
survey for about 1/4 sky area. The SDSS data have been 
made public in a series of yearly data releases and the 
present is the Seventh Data Release (DR7), including data 
up to the end of SDSS-II in July, 2008. 

The optical spectrum of RE J1034-I-396 is downloaded 
from SDSS DR7. Its spectrum covers the wavelength range 
of 3800-9200 A with a spectral resolution of 1800 < R < 
2100. At its redshift of 0.043, the projected fiber aperture 
diameter of 3" is about 2.4 kpc, containing most light from 
its bulge. With nuclei spectrum superimposed on it, the stel- 
lar absorption features in its SDSS spectrum provide us the 
possibility to investigate the property of its host bulge, in- 
cluding the bulge stellar velocity dispersion, the stellar mass, 
and the star formation history. We did not apply aperture 
correction to the stellar velocity dispersion because this ef- 
fect can be omitted for z < 0.3 (Bernardi et al. 2003). 

We outline our steps to do the SDSS spectral analysis. 
(1) We use SSP synthesis (STARLIGHT; Cid Fernandes et 
al. 2005) to model the stellar contribution in the Galactic 
extinction-corrected spectrum in the rest frame (Cid Fernan- 
des et al. 2005; Bian et al. 2006, 2007, 2008). The Galactic 
extinction law of Cardelli, Clayton & Mathis (1989) with 
Rv = 3.1 is adopted, which is also used for the host ex- 
tinction with V-band extinction (from mag to 5.0 mag). 
We use 45 default templates in Cid Fernandes et al. (2005), 
which are calculated from the model of Bruzual & Char- 
lot (2003). The linear combination of 45 templates is used 
to represent the host bulge spectrum. These 45 templates 
comprise 15 ages, t = 0.001, 0.00316, 0.00501, 0.01, 0.02512, 
0.04, 0.10152, 0.28612, 0.64054, 0.90479, 1.434, 2.5, 5, 11 
and 13 Gyr, and three metallicities, Z — 0.2, 1 and 2.5 Zq 
(Cid Fernandes et al. 2005). At the same time as the SSP 



fit, we add a power-law component in the code to represent 
the AGN continuum emission, and an optical Fe ll template 
from the prototype NLSl I ZW 1 (Boroson & Green 1992) 
to model the Fe ll emission. We exclude the AGN mission 
lines, such as H Balmer lines, [O ll] A3727, [Ne ill] A3869, [O 
III] AA4959, 5007, [N ii] AA6548, 6583, [S ii] AA6717, 6731. 

The synthetic spectrum is built using the following 
equation. 



E 

.j=i 



(1) 



where bj^x is the j^^ template normalized at Ao = 4020 A, Xj 
is the fiux fraction at 4020 A, Mxq is the synthetic fiux 4020 
A, rx = 10"° ''''*^ --^Ao) ig the reddening term by V-band ex- 
tinction Av, and G{vo,Vd) is the line-of-sight stellar velocity 
distribution, modeled as a Gaussian centered at velocity no 
and broadened by the velocity dispersion Vd- Due to differ- 
ent velocity dispersion in the stellar lines and Fe ll lines, we 
use another line-of-sight Fe ll velocity distribution G{vfe, Ufe) 
for the Fe ll emission. X[a is the Fe ll fiux-fraction at 4020A. 
The line-of-sight Fe ll velocity distribution G{vic, f^fc) is also 
modeled as a Gaussian centered at velocity vtc and broad- 
ened by the velocity dispersion afo. We first fit the Fe ll lines 
and the continuum in the fitting windows. And we find that 
ate = 411 ± 46km s~^. ate is fixed by 411 km/s in the fitting 
of SSP and Fe ll . When we change crfo, the SSP and Fe ll 
fitting results do not change, considering the errors. We also 
exclude the Fe ll emissions and do the SSP fit, considering 
the errors, the SSP results do not change. The best fit is 
reached by minimizing reduced x^, 



x\x, Mxo,Av, vo,vd)=Y^ [{Ox - Mx) wxf /N 



(2) 



where the weighted spectrum wx is defined as the noise as- 
sociated with each spectral bin as reported by the SDSS 
pipeline output, N is the total unmasked pixels, is calcu- 
lated by the difference between observed spectrum and the 
model spectrum in the fitting of SSP and Fe ll . For RE 
J1034-H396, the best fit of SSP and Fe ll gives = 1-4 (Fig 
1). We find that Ay = 0, and the host extinction can be 
neglected (Fig 1). Through above spectral synthesis, we can 
obtain some parameters, such as bulge velocity dispersion 
Vd, fiux-fraction Xj, mass- fraction pj, stellar mass i\f, . The 
mass-flux ratio can be found in STARLIGHT manual on 
the web site ^http: / /www. sta rlight. ufsc.brT] The results are 
shown in Fig. 1 and Table 1. 

With the simulation, it is suggested that the uncertainty 
of SSP results can be given by the effective starlight signal- 
to- noise (S/N) at 4020 A (Cid Fernandes et al. 2005; Bian 
et al. 2007). The S/N at 4020 A is 37 and the starlight 
fraction at 4020 A is about 43% (Table 1). Therefore the 
effective starlight S/N at 4020 A is about 16, corresponding 
to an uncertainty of 8 km s~^ for the velocity dispersion, 7% 
for the mass-fraction, 8% for the flux-fraction, 0.06 dex for 
stellar mass logMt (Cid Fernandes et al. 2005; Bian et al. 
2007). We adopt 10% as the uncertainty of host bulge flux 
fraction. 

(2) Considering the broad wing in the H/3 line profile, 
two broad components are used to model broad H/3 pro- 
file from BLRs. The H/3 line from the narrow line regions 
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Figure 1. Results of simple stellar population (SSP) synthesis for the SDSS spectrum of RE J1034+396. Top panel: the left panels 
are the observed and the synthesis spectra (black and red, respectively) , and the residual spectrum. The panels on the right show the 
the flux-fraction Xj and mass-fraction fij versus the ages of the simple stellar populations. Bottom panel: the black line is the observed 
SDSS spectrum corrected for Galactic extinction in the rest frame. The cyan line is the power-law continuum. The green line is the host 
contribution. The blue line is the Fe ll spectrum. The magenta line is the residue. The red line is the composition of the host contribution, 
Fe II lines and the power-law continuum. 
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Table 1. The host bulge properties of RE J 1034+396 



(T, (km s ^) 


logM,{MQ) 




/host 


(1) 


(2) 


(3) 


(4) 


67.7 ±8 


10.155 ±0.06 


0.0163 ±0.0011 


43 ± 4% 



Note. Col(l); the host bulge stellar velocity dispersion; Col(2): 
the stellar mass; Col(3): the host bulge specific star formation 
rate; Col(4): the host bulge flux fraction at 4020 A with 10% 
uncertainty. 



Table 2. The emission lines properties of RE J1034+396 





FWHM" 


flux* 


H/3 


1690 ± 296 


1.1 ±0.84 




617 ± 155 


0.88 ±0.73 




1022 ± 4592" 


0.10 ±0.63 




292 ± 114" 


0.12 ±0.13 


[O III] 


1023 ± 18" 


1.26 ±0.03 




292 ± 5" 


1.43 ±0.03 


Ha 


969 ± 552 


3.98 ±3.15 




3804 ± 1100 


2.25 ±0.39 




272 ± 25" 


1.16 ±0.25 




1162 ± 275" 


1.09 ±2.75 


[N ii] A6583 


278 ± 80" 


1.16 ±0.08 




1162 ± 287" 


1.1 ±2.81 


Fe II 


967 ± 108 


1.76 ±0.2= 



Note. °: FWHM of a Gaussian proflle in units of km s~^ . in 
units of 10^^* erg s~^ cm~^. Fe II flux between 4434 A and 
4684 A. ": the components from NLRs. The second weaker broad 
component of Ha is probably due to the contribution from the 
continuum subtraction. 



(NLRs) has the same profile as the [O ill] line from NLRs. 
Two components are used to model the asymmetric [O ill] 
line profile, and two components are used to model H/3 line 
profile from NLRs (Bian et al. 2008). Therefore, total four 
Gaussians are used to model the H/3 line profile, and two sets 
of two Gaussians are used to model the [O ill] AA4959, 5007 
lines. As the same as the H/3 profile, four Gaussians are used 
to model the Ha profile (two broad components from BLRs 
and two narrow components from NLRs), and two sets of two 
Gaussian are used to model the [N ll] AA6548, 6583 lines. We 
take the same line width for each corresponding component 
of [O III] AA4959, 5007 and H/3 from NLRs, fix the flux ra- 
tio of [O III] A4959 to [O lll] A5007 to be 1:3, and set the 
wavelength separation to the laboratory value. We take the 
same line width for each corresponding component of [N ll] 
AA6548, 6583 and Ha from BLRs, fix the flux ratio of [N ll] 
A6548 to [N ll] A6583 to be 1:3, and set the wavelength sepa- 
ration to the laboratory value. For RE J1034+396, the best 
fits of the H/3 , Ha lines give reduced — 2.2, = 2.14, re- 
spectively. The results are shown in Fig. 2 and Table 2. The 
third component of the H/3 line from NLRs is weaker, the 
results do not change if we remove this weaker component. 




6450 6500 6550 6600 



Rest Wavelength (A) 

Figure 2. An example of the line fit for RE J1034+396. Top 
panel: the line fit for H/3 and [O III] line. The black line is the orig- 
inal spectrum after Galactic-extinction correction, the starlight 
subtraction, the power-law continuum and the Fe II subtraction 
in the rest frame. The residual is shown at the bottom. The mul- 
tiple Gaussian components are in blue and the sum of them is in 
red. The fitting window is in green. Bottom panel: the line fit for 
Ha and [N ll] line. The line colors are the same to that in top 
panel. 

3 RESULT AND DISCUSSION 
3.1 Mass 

Firstly, we derive the mass from the bulge stellar velocity 
dispersion, a,. By SSP method, the measured bulge stellar 
velocity dispersion Vd is 28 ± 8 km s~^ . Considering the 
resolutions of the SDSS spectra and the template spectra, 
we have to do these corrections by the following formula Q 

^* = \J'"l+ '^temp - ^Lst (3) 

where the SDSS spectral resolution (Tinst is adopted as 60 
km s~^ (Greene & Ho 2005a), the template spectral res- 
olution (Jtomp is adopted as 86 km s~^ (Cid Fernandes et 
al. 2005). Assuming that the errors on o"inst and o"temp are 
negligible, the error on cr, will be the same as that on Vd- 
iJd = 28 ± 8km s"^ leads to cr* = 67.7 ± 8km s"^ Using 
the Mbh — cr* relation for quiescent galaxies, Mbh(o'*) = 
10** "[(j./(200 kms~^)]* °^ Mq (Tremaine et al. 2002), we 

^ |http://www.starlight.ufsc.br/| 
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derive the SMBH mass is (1.7+,!,;?) x 10** Mq for cr, = 
67.7 ± 8km s~^. The uncertainty of 8 km s^^ for a, would 
lead to an error of 0.1 dex for log Afan . The total uncer- 
tainty is about 0.32 dex considering the error of 0.3 dex 
from the Mbh — cr* relation (Tremaine et al. 2002). There- 
fore, the SMBH mass of RE J 1034-1-396 derived from cr, is 
about (1 - 4) X 10*^ Mq. 

Secondly, we derive the mass from the broad H/3 second 
moment cth/j and the empirical _Rblr — ALa(5100A) relation 
(Kaspi et al. 2000, 2005; Bentz et al. 2006; Bian et al. 2008). 
With the power-law component, we can derive the nuclear 
flux at 5100 A is 45 X 10"" erg s"^ cm-^ and ALa(5100A) 
is 5.6 X lO'*^ erg s~^. The first moment of the line profile 



P(A) is 



Ao = 



/ AP(A)dA 
/ P(A)dA ' 



The second moment of the line profile P(A) is 



2 



JX^P(X)dX 

J P{\)d\ 



(4) 



(5) 



For a Gaussian line profile, the ratio of FWHM to the second 
momentum is FWHM/criinc = V81n2 2.35; while for a 
Lorentzian profile, aiinc ^ oo. For RE J1034-I-396, the H/9 
second momentum cth/j is calculated from the reconstructed 
H/9 profile of two broad H/3 components from BLRs in step 
2 of spectral analysis, and cth/s is 577 ± 144 km s~^ . From 
the reconstructed H/3 profile from BLRs, FWHM is 802 ± 
200km s"\ The SMBH mass can be calculated from the H/? 
second moment by the following equation(Bian et al. 2008) , 



Mbh = / 



PblrAF" 



G 



f X 7.629 X 10' 



6r AL;,(5100A) o.5i8r 

lO^-ierg s-i ^ 



0"H/3 



lOOOkm s 



where Rblr = 39.08 x [ALA(5100A)/(10**erg s-i)]°-5i** 
hght-days (Bentz et al. 2006), and / = 3.85 (Colhn et al. 
2006). With cTH/s = 577 ± 144km s'^ and / = 3.85, the esti- 
mated SMBH mass is (0.7 - 3.4) x lO'^M©. The uncertainty 
of the mass calculation from the H/3 line is mainly from 
the systematic uncertainties, up to about 0.5 dex, which is 
due to the unknown kinematics and geometry in BLRs (e.g. 
Krolik 2001; Peterson et al. 2004). Grupe et al. (2004) gave 
ALa(5100A) as 1.5 x 10*^ erg s"^ Considering the half light 
contribution from the host bulge, the nuclear XL\{5100A) 
is consistent with ours. They also gave the H/? FWHM as 
700 ± 110 km s"^ . When using this FWHM value to cal- 
culate the mass, the mass would be smaller by 0.5 dex than 
that from the cth/j (Bian et al. 2004, 2008). The new mass 
from the H/3 second moment is consistent with that from the 
Mbh — o"* relation. Hereafter, we adopt (1 — 4) x 10^ Mq as 
its SMBH mass. This SMBH mass value of RE J1034-I-396 
is consistent with that derived from the accretion disk fit- 
ting of its SED (Puchnarewicz et al., 2001; Wang & Netzer 
2003). 

It is suggested that the gaseous kinematics of NLRs 
be primarily governed by the bulge gravitational potential. 
Using the linewidth of [O in] or [N ll] to trace the a,, we 
find that, from narrow [O ill] component, cr[oiii] = 124 km 
s~ (see Table 2). It would lead to a mass of 2 X IO'^Mq, 
which is larger than mass from a, and (th/3 by an order of 



magnitude. For RE J1034-I-396, crjoin] cannot be used to 
trace a*. 



3.2 High-frequency QPOs and BH spin 

QPOs can be classified as low-frequency QPOs (roughly 
0.1-30 Hz) and high-frequency QPOs (roughly 40-450 Hz, 
Rimillard & McClintock 2006). Assuming QPO frequency 
(j^bhb) scales inversely with the BH mass (Mbhb) in BHBs, 
A/bh/Mbhb = vbub/vbh, where vbh is the QPO frequency 
for SMBH in AGN. By vbhb = I'BH X {Mbh/Mbhb), the 
QPO frequency of 2.7 x 10"* Hz in RE J1034-I-396 (its 
SMBH mass of 2 x 10^ M© with an uncertain of a factor 
2) corresponds a frequency of 27-108 Hz in Galactic BHBs 
with a 10 Mq BH. Our resuh of 27-108 Hz suggests that the 
QPO found in RE J1034-I-396 belongs to a high-frequency 
QPO. Larger estimated SMBH mass, smaller adopted BHBs 
mass would make this frequency value larger. 

It is believed that the high-frequency QPOs may depend 
only on the mass and the spin of the BH. This dependence 
can be expected for coordinate frequencies (see Merloni et 
al. 1999) or for disk oscillation modes in the inner accretion 
disk (see Kato 2001). For a Schwarzchild BH, the inner- 
most stable circular orbit (ISCO) corresponds a maximum 
orbit frequency as uisco ~ 22OO(A/bh/A/0)"^Hz. For an ex- 
treme Kerr BH, i^sco = 1615O(Mbh/M0)-^Hz (Rimillard 
& McClintock 2006). For three BHBs with high-frequency 
QPOs and well-constrained BH masses, Rimillard & Mc- 
Clintock (2006) suggested an empirical frequency - mass re- 
lation for the high-frequency QPOs, vo — 931(Mbh/M0)~^ 
Hz. With the SMBH mass of 2 x 10*^ Mq in RE J1034H-396, 
the corresponding frequencies derived from the ISCO of a 
Schwarzchild BH, the ISCO of an extreme Kerr BH, and 
2 the empirical frequency - mass relation are 11.0 x 10~*, 

-] Mq (6) g Q ^ ^Q-i^ 4 7 ^ ^Q-i jj^^ respectively. Considering an un- 



certainty of a factor of 2 in the SMBH mass and other uncer- 
tainty in empirical frequency - mass relation for the high- 
frequency QPOs (e.g., the BH spin), the double or triple 
value of the fundamental frequency 4.7 x 10"* Hz dose not 
deviate much from the observed frequency of 2.7x 10"'' Hz in 
RE J1034-I-396. It implies that the QPOs phenomenon has 
the same origin in BHBs and AGN (e.g., Gierlinski 2008). 



3.3 Star formation 

The widely employed method to measure the current star 
formation rates (SFR) is by means of the ultraviolet con- 
tinuum, the emission lines of Ha , [O ll] , et al. (Kennicutt 
1998). For a sample of 82302 star formation galaxies from 
SDSS, Asari et al. (2007) investigated the current SFR de- 
rived from the SSPs synthesis and found that it is consistent 
with that from the Ha SFR indicator (see their Fig. 6). The 
current specific SFR (SSFR) is defined by the mean SSFR 
over the past O.lGyr, i.e., first 21 of the 45 templates de- 
scribed in section 2 with age less then 0.1 Gyr (Asari et al. 
2007; Chen et al. 2009), 



SSFR(t < O.lGyr) = 



E21 
_ j^i W 

A-f. ^ O.lGyr 



(7) 



For RE J1034-K396, current SSFR(t < O.lGyr) is 0.0163 ± 
0.0011 Gyr~\ logAf. is 10.155 ± 0.06Afo , and current SFR 
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is 0.23±0.016 MQyr-^. Using the SFR formula from Kenni- 
cutt et al. (1998), this current SFR of 0.23 ± 0.016 M^yr'^ 
corresponds to the Ha luminosity as (2.9±0.2) x 10'*° erg s~^ 
or the [O ll] luminosity as (1.6 ±0.1) x 10'"^ erg s~^ coming 
from the starburst process. From the SDSS spectrum, we 
find that total Ho? luminosity and total [O ll] luminosity is 
4.4 X 10*^ erg s~^ , 3.4 x 10*° erg s"^ . Therefore, about 
6.1%-7.0% of the total Hq luminosity and 44%-50% of the 
total [O ll] luminosity come from the star formation (Ho 
2005). 

We also use 150 templates instead of 45 templates in 
SSP synthesis to model the stelar absorption contribution, 
like in Asari et al. (2007). We found that a* = 72.4km s~\ 
logM. = 10.244, SSFR(t < O.lGyr) = 0.02786Gyr-\ Using 
150 templates, the cr, and logM, are consistent with that 
using 45 templates, the stellar flux fraction is 47%, and the 
SSFR would be larger by 70%. 

Using type-H AGN sample from MPA/JHU catalog 
(Kauffmann et al. 2003), Chen et al. (2009) found a cor- 
relation between the Eddington ratio A (the ratio of the 
bolometric luminosity, Lboi, to the Eddington luminosity, 
I/Edd) and the mean SSFR, logA = (-0.73 ± 0.01) + (1.5 ± 
0.01)log(SSFR/Gyr^^), suggesting that supernova explo- 
sions play a role in the transportation of gas to galactic cen- 
ters. Considering that the difference between type H AGN 
and type I AGN is due to the orientation of the line of sight, 
type I AGN would also follow this logA — logSSFR relation. 
Here we calculate the Eddington ratio for RE J1034-I-396. 
Using Lboi = 9 X ALa(5100A) (Kaspi et al. 2000), I/Edd = 
1.26 X lO^*(Mbh/M0) erg s^\ and Afen = 2 x 10® Mq, we 
find that the Eddington ratio is 0.2. Considering the un- 
certainty of a factor of 2 in A/bh , the uncertainty of the 
Eddington ratio is a factor of 2, i.e., logA = — 1 ~ —0.4. 
Using the formula of (3) in Chen et al. (2009), the SSFR of 
0.0163 ±0.001 Gyr"^ leads to logA as -3.45 ~ -3.37, which 
deviates much from the value of — 1 ~ —0.4. The deviation 
of logA is 2.37 ~ 3.05 in the logA — logSSFR diagram of Chen 
et al. (2009). Considering most of the bolometric luminos- 
ity is emitted in the UV and soft X-ray band, Grupe et al. 
(2004) estimated the bolometric luminosity from a combined 
power-law model fit with exponential cutoff to the optical- 
UV data and a power law to the soft X-ray data. They found 
Lboi ~ 22 X ALa(5100A), which is consistent with its hot big 
blue bump. If that is the case, logA would be larger by 0.39. 
Therefore, RE J103 4±396 does not follow the relation be- 
tween the Eddington ratio and SSFR found by Chen et al. 
(2009), although a large scatter in their relation. The reason 
for the deviation would be explored in other place. 



4 CONCLUSIONS 

The host bulge of a NLSl RE J1034±396 with X-ray quasi- 
periodic oscillations is investigated through its optical spec- 
trum from SDSS DR7. The main conclusions can be sum- 
marized as follows: (1) The host bulge flux contribution at 
4020 A is about (43 ± 4)% in the SDSS spectrum and the 
bulge stellar velocity dispersion cr, is 67.7 ± 8 km s~^ . 
Considering the scatter in A^bh — o", relation, the SMBH 
mass from AfBH - cr, relation is (1 - 4) x lO'^Af©. (2) Us- 
ing multi-Gaussians to model the H/3 , [O ill] , Ha , and 
[N ll] emission lines, we find that the mass from cth/s and 



ALa(5100A) is (0.7 ±3.4) x 10*^ Mq, consistent with that 
from a*. However, the [O ill] velocity dispersion is 124±2 km 
s^^ , about a double of cr, . (3) The SMBH mass of about 
(1 - 4) X IO^A/q implies that the QPO of 2.7 x IQ^^Hz 
found in RE J1034±396 can be scaled to a high-frequency 
QPO at about 27-108 Hz found in Galactic BHBs with a 10 
Mq BH. (4) From the SSPs synthesis, we find that current 
SSFR(t < O.lGyr) is 0.0163 ± 0.0011Gyr-\ \og{M^/MQ) 
is 10.155 ± 0.06, and current SFR is 0.23 ± 0.016 A^0yr"^ 
RE J1034+396 does not follow the relation between the Ed- 
dington ratio and SSFR. About 6.1%-7.0% of the total Ha 
luminosity and 44%-50% of the total [O ll] luminosity come 
from the star formation process. 
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